Dynamic first-order phase transition in kinetically constrained models of glasses 
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We show that the dynamics of kinetically constrained models of glass formers takes place at a first- 
order coexistence line between active and inactive dynamical phases. We prove this by computing 
the large-deviation functions of suitable space-time observables, such as the number of configuration 
changes in a trajectory. We present analytic results for dynamic facilitated models in a mean-field 
approximation, and numerical results for the Fredrickson- Andersen model, the East model, and 
constrained lattice gases, in various dimensions. This dynamical first-order transition is generic in 
kinetically constrained models, and we expect it to be present in systems with fully jammed states. 

PACS numbers: 05.40.-a, 64.70.Pf 
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An increasingly accepted view is that the phenomenol- 
ogy associated with the glass transition [l| requires a 
purely dynamic analysis, and does not arise from an un- 
derlying static transition (see however Q). Indeed, it has 
been suggested that the glass transition manifests a first- 
order phase transition in space and time between active 
and inactive phases |3|. Here we apply Ruelle's thermo- 
dynamic formalism [J, Q to show that this suggestion 
is indeed correct, for a specific class of stochastic mod- 
els. The existence of active and inactive regions of space- 
time, separated by sharp interfaces, is dynamic hetero- 
geneity, a central feature of glass forming systems @|. 
This phenomenon, in which the dynamics becomes in- 
creasingly spatially correlated at low temperatures, arises 
naturally p} in models based on the idea of dynamic facil- 
itation, such as sp in-facilitated models [1, Q , constrained 
lattice gases flOl . Fill] and other kinetically constrained 
models (KCMs) pjj • Fig. Q] illustrates the discontinu- 
ities in space-time order parameters at the dynamical 
transition in one such model, together with the singular- 
ity in a space-time free energy, as a function of a control 
parameter to be discussed shortly. 

The thermodynamic formalism of Ruclle and cowork- 
ers was developed in the context of deterministic dynam- 
ical systems Q. While traditional thermodynamics is 
used to study fluctuations associated with configurations 
of a system, Ruelle's formalism yields information about 
its trajectories (or histories). The formalism relies on the 
construction of a dynamical partition function, analogous 
to the canonical partition function of thermodynamics. 
The energy of the system is replaced by the dynamical 
action (the negative of the logarithm of the probability 
of a given history); the entropy of the system by the 
Kolmogorov-Sinai entropy [l3j . and the temperature by 
an intrinsic field conjugate to the action. This formal- 
ism has been exploited recently to describe the chaotic 
properties of continuous-time Markov processes [f| • 

In this work, we define the dynamical partition sum 0, 
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FIG. 1: First order transition in terms of the field s. (Top) 
The dynamical order parameter K(s) (the average number of 
configuration changes in a trajectory) and its large-deviation 
function ^k(s) for the FA model, calculated in a mean-field 
approximation, for d = 3 and T — 0.5: see Eqs. (|6I9|I . The 
large deviation function is singular at s = and the order 
parameter K has a first-order jump. The dynamics has two 
phases, an active one for s < and an inactive one for s > 0. 
Physical dynamics take place at s = 0, where the two dy- 
namic phases coexist. (Bottom) An alternative order param- 
eter pk{s) (the average number of excited sites in a trajectory: 
see Fig. [2} in the d — 1 FA model at T = 0.91, calculated nu- 
merically in a finite system (N = 100 sites). The transition 
is absent when the kinetic constraints are removed. 



IH for our stochastic systems by 

Zx{s,t) = Prob(history)e 



— sX(history) 



(1) 



histories 



where the sum is over histories from time to time t; the 
probability of a history is Prob(history); and ^"(history) 
is the number of configuration changes in that history. 
iC(history) is a direct measure of the activity in a history: 
an active trajectory has many changes of configuration, 
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FIG. 2: Large-deviation function %pK and order parameters as a function of system size, for the FA model (d = 1), the 
TLG (d = 2) and the East model (d = 3). For the FA and East models, we show the average density of excited sites 
p K {s) = Z( S ,t)- 1 E histori c S Prob ( histoi 'y)e" sj5r(history) p(history), where p(history) = (iVi)" 1 /Jdt' m(t'). For the TLG, we 
replace p(history) by r(history) = (Nt)" 1 f Q dt' r(C(t')), to obtain the averaged escape rate r_rc(s). We use representative 
conditions: T = 0.91 for the FA and East models, and density 0.5 in the TLG. Results were obtained by the dynamical method 
of[IJ. 



an inactive one has few or none. At long times 

Z K (s,t) <x e^ s \ (2) 

where 4>k(s) is the large deviation function for K. The 
quantity —J^k(s) is a free energy per unit time, for tra- 
jectories [lj]. When t is large, the derivatives of iPk(s) 
give the cumulants of K. For example, the mean num- 
ber of configuration changes in a trajectory of length t is 
A"(0), where 

K(s) = -t^ K (s). (3) 
as 

For some models, it is reasonable to use the dy- 
namical action, or dynamical complexity, Q + = 
In [Prob (history)] as an order parameter for dynami- 
cal activity, in place of K. In that case the large devia- 
tion function is the topological pressure and the average 
of (—Q+/t) at s = is the Kolmogorov-Sinai entropy 
d, |||. In 0|, the action was denoted by £ = —Q+ and 
the quantity 1 + s was denoted by b. In the following, we 
focus on K, which simplifies the analysis. 

We consider facilitated models |12| defined by a bi- 
nary field rii = 0, 1 and a non-interacting Hamiltonian 
H = J^i n ii where i = 1, . . . , N are the sites of a lat- 
tice. The dynamics at site i is subject to a kinetic con- 
straint Ci{rij) which is a function of the neighbours rij 
of i. That is, site i changes its state with a rate propor- 
tional to Ci. For example, in the Fredrickson- Andersen 
(FA) model [1, [HJ , Ci = 1 if any of the nearest neigh- 
bours j of i are in the "excited" state, rij = 1; otherwise 
Ci = 0. Sites with Cj = 1 make the transitions — ► 1 
with rate c and 1 — > with rate (1 — c), where c is 
the equilibrium concentration of excitations at tempera- 
ture T, c = (rii) = (1 + e 1 / T )~ 1 . One can also consider 



exchange dynamics l^Oj <-* Ojlj between nearest neigh- 
bour sites i and j, so that the total density of excited 
sites is conserved. The result is a constrained lattice gas 
[ToL [Til [l2j • For example, the two vacancy facilitated lat- 
tice gas, or (2)-TLG [ll(, is defined on a triangular lat- 
tice, with exchange dynamics and a constraint dj which 
vanishes unless the two common nearest neighbours of 
sites i and j are vacant. 

Dynamic heterogeneity in KCMs is a consequence of 
space-time correlations in the dynamics. Within the 
glassy regime, trajectories contain space-time regions 
with relatively fast dynamics, and others which rearrange 
slowly 7]. Slow regions can be non-finite in extent, that 
is, KCMs contain states for which the number of sites 
with Ci ^ is subextensive in the system size N. These 
configurations are important for our purposes. Transi- 
tions out of these states occur with a rate that is subex- 
tensive in N. For systems obeying detailed balance, we 
show below that if any such states exist the partition 
sum ([TJ is dominated by trajectories localised in them, 
at large N and t and for all positive s. The order pa- 
rameter K is subextensive for these trajectories, so that 
K(s)/N vanishes for s > in the large N limit. Con- 
versely, for s < 0, the dominant trajectories visit states in 
which a finite density of sites have Ci ^ 0, and K(s)/N 
is finite. These arguments establish the existence of the 
transition shown in Fig. [TJ 

For a more quantitative analysis, we identify ^k(s) as 
the largest eigenvalue of a time evolution operator 0, [l5| . 
The dynamics of the stochastic model are specified by a 
master equation 

d t P(C, t) = -r(C)P(C, t) + J2 W(C -> C)P(C, t), (4) 

a 

where C denotes a configuration of the system, the 
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W(C — > C) are the transition probabilities between con- 
figurations, and r(C) — J^ c , W(C — ► C) is the rate of 
escape from C. We write this equation in an operator 
form dtP = WP, where the matrix elements of W are 
W(C,C') = W{C ->C)-r(C)5cfi>. 

Now, consider the probability P(C,t,K) that the 
system is in configuration C at time t, having made 
K transitions since time 0. We define P(C,t,s) = 
^e^fC,^), so that Z(s,t) = £ c P(C,*,s)- The 
time evolution of P obeys a Master equation d t P — 
W K P, where 0,[H 



<K 



{C,C) = e~ s W{C ^ C) ~ r{C)5 C fi> 



(5) 



The calculation of tPk(s) reduces to finding the largest 
eigenvalue of Wjj. 

The FA model dynamics can be described using the 
Doi-Peliti [16j occupation number formalism. The oper- 
ator ||SJ) then reads [13, El: 



4 FA) 



E'< 



[e(e 



(6) 



where aj and a,i are bosonic creation and annihilation 
operators, hi — a\a,i and indicates that the sum is 
over pairs of nearest neighbours. For simplicity we have 
allowed for multiple occupancy per site: changing the 
dynamics in this way has little effect at low temperatures 
hj. The single occupancy case can be formulated in a 
similar manner [l8j . and the results are analogous. 



We begin with a simple mean- field treatment of ^ 



4 FA) - 



We discard fluctuations and define wj^ A ^ ((f), <fi) by mak- 



ing the replacements a 
in the operator W^ A ^ . 

4 FA) 



eigenvalue of 



; — > (a*) = (f) and aj — > (aj) = <fi 
Then we can estimate the largest 
by finding the stationary points of 



W 



, (f>) . The Euler-Lagrange equations are 

= 2# (<f) ~ e- s c) + <j> (c - e- s 4>) , (7) 
= 2(0-e- s )0(/. + c0(l-e- s 0) . (8) 



There are two solutions that correspond to steady states 
of the dynamics: (i) <j) = aj) = c<f)^ with 0„ = |e _s 



w 



8; and (ii) (f> = <f) = 0. For solution (i), 
is positive for s < and negative otherwise; 



( FA )^ 1\ 



0. Our estimate for the 



for solution (ii), 
large deviation function is given by the larger value of 

Wjr A) (&0), so that 

V,^- f -)( s ) = / Nd ( c ^) 2 ( e_ ^ - !) ( s < °) (9) 
Vk K ' \ (s > 0) V ; 

where d is the spatial dimension, and 0„ is defined above. 

Figure DJtop) shows tpK ^ V s ) anc ^ K( s )- The active 
phase is stable for s < 0; the inactive one is stable for 
s > 0. The large deviation function V'k^M ^ s con t mu_ 
ous at s = 0, but its derivative K(s) displays a first-order 
jump. 
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FIG. 3: (a) Space-time phase diagram for facilitated models. 
s = is a first-order transition line between the active phase 
(s < 0) and the inactive phase (s > 0). It terminates in a 
critical point at T = 0. (b) Mean-field variational free-energy 
Tk{p,s), Eq. |TT}, for the FA model at T = 0.5, for several 
values of s. For s < the active phase is dominant and 
the absolute minimum of Tk is at p 7^ 0. At s — there 
is dynamic phase-coexistence. For s > the inactive phase 
dominates and p — minimizes the free-energy. The active 
minimum no longer exists beyond the spinodal s sp . 



The mean-field approximation is clearly very crude. 
It can be improved systematically with loop corrections, 
but the coexistence scenario does not change qualita- 
tively: it is a consequence of the kinetic constraints and 
the existence of states with subextensive escape rates. 
We now confirm this by means of numerical simulations. 

The time evolution associated with the operator Wk 
can be obtained dynamically. This operator does not 
conserve probability [see Eq. (JSj)] but there are methods 
to simulate the evolution that it represents p^ |. In our 
case we apply this scheme, analogous to the quantum dif- 
fusion Monte Carlo procedure, to continuous time Monte 
Carlo dynamics [2fJ. In Fig. [1] we show a plot of pk{s) 
for the FA model in one dimension. In Fig. [5] we present 
a finite size scaling analysis of both i^k and a dynamical 
order parameter for the FA model in one dimension, the 
(2)-TLG model in two dimensions, and the East model in 
three dimensions (see e.g. [HI). As we approach N — > 00 
the crossover between active and inactive phases becomes 
sharp. Figure shows results for three different models, 
in different dimensions and with different dynamical con- 
straints; the TLG has a conserved density while the FA 
and East models do not. Our similar results for these 
different models demonstrate that space-time phase co- 
existence occurs quite generally in KCMs. 

We have demonstrated that, for any temperature T, 
the dynamics of KCMs such as the FA or East models 
take place at dynamical phase coexistence between active 
and inactive phases. We can construct a phase diagram 
in terms of T and s, see Fig. [3][a). The s = axis is a 
first-order transition line. It ends in a critical point at 
T = 0. For constrained lattice gases such as the Kob- 
Andersen model and the (2)-TLG, the phase diagram is 
similar, with T replaced by density of particles c. In that 
case the critical point is at the maximum density c = 1. 
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We emphasize that while a zero temperature dynami- 
cal critical point is common to many KCMs, this is not a 
sufficient condition for space-time phase coexistence. For 
example, consider the pair appearence and annihilation 
(AA) model of [18(, which has the same critical proper- 
ties as the FA model. Its dynamical rules are %A ^ A$, 
AA — > 0, 00 — > AA, with rates D, A and 7, respectively. 
All states in this model have extensive escape rates, so 
we do not expect any transition at s = 0. Deriving the 
Euler-Lagrange equations analogous to (|7I8[) . leads to a 
large deviation function that is analytic at s = 0. We 
have also calculated tPk{s) exactly [22| for the AA model 
in d = 1, at the free fermion point A + 7 = 2D [23| : 
the large deviation function ipn(s) is indeed analytic at 
s = 

Finally, in models that obey detailed balance with re- 
spect to a probability distribution pc, such as the ones 
considered here, iPk{s) can be calculated through a vari- 
ational method. The master operator Eq. |(5J) is made 
symmetric by a similarity transformation, Mk(C,C) = 

Pc 1/2 W jr (C,C')Pc' 2 = e- 3 y/W(C -+ C)W(C C) - 
r(C)6c.c> ■ Since Hr- is symmetric and has the same eigen- 
values as Wjf, we can apply a variational principle: 



iI>k(s) = max 

Vc 



j: c VcM K (C,C')Vc 
Ec We 



(10) 



Using a trial distribution in which only one of the Vc 
is finite shows that the largest diagonal element of H is 
a lower bound on ^r-(s), so that 4>k{s) > — mine r(C). 
Hence, if there exists a state for which the escape rate 
r(C) is subextensive, then ipK(s)/N > at large N. Fur- 
ther, iPk{s) is non-increasing [since K is non-negative] 
and i/jk(0) = 0, so the existence of a subextensive escape 
rate establishes immediately that iPk{s)/N — K(s)/N = 
for all s > in the large N limit, as we asserted above. 



As an example of how to obtain iPk{s) from Eq. (fT0|) 
consider the FA model in a mean-field geometry, such as 
the complete graph. The transition rates in this case are 
W(n n + 1) = cn, Win + 1 n) = n(n + 1)/N, 
where n = rii is the total number of excitations. For 
the variational state we assume V n = e N ^ n / N \ for some 
function f(p) of the excitation density p = n/N. In the 
limit of large N, the leading contribution to (fT0|) comes 
from values of p that maximize f(p), and Eq. (|10[) reduces 
to ipi<{s) = — min p J-k(p, s), where 



F K (p,s)=Np (p -2e 



cp 



(11) 



Tk is a Landau free-energy function for the order param- 
eter p. The minimum of Tk occurs at p* = c(0 + ) 2 , and 
from Eq. (fTU|) the result of Eq. ^ is recovered. As shown 
in Fig. [3th) the function Tk(p, s) behaves in the charac- 
teristic way associated with a first-order phase-transition. 

We have shown that the dynamics of KCMs is char- 
acterised by the coexistence in space-time of active and 
inactive dynamical phases. In our view, this dynami- 
cal phase coexistence underlies the heterogeneous parti- 
cle dynamics observed in glass formers. Thus, experi- 
mentally obervable phenomena such as transport decou- 
pling pi] arise from the fluctuations associated with this 
dynamic phase equilibrium. The extension Q of our re- 
sults to atomistic models will clarify the degree to which 
the theoretical framework described here captures the 
generic features of glassy dynamics. 
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